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ABSTRACT 


This paper presents a thorough exploration of vehicle dynamics modeling and simulation, focusing on lateral, longitudinal, and 
vertical motion. A comprehensive mathematical representation of a car-like vehicle is developed, incorporating both kinematic and 
dynamic models to accurately capture its behavior. Initially, a kinematic model is established to describe fundamental motion in a 
2D plane, while a subsequent dynamic model considers various forces influencing the vehicle's motion in a 3D space. The vertical 
motion of the vehicle is primarily influenced by its suspension systems. In this paper, we initially introduced a fundamental 2 Degree 
of Freedom (DOF) quarter car suspension model and subsequently presented a more comprehensive half car suspension model 
which accounts for both the front and rear body parts of the vehicle. The developed models are implemented using 
MATLAB/Simulink, enabling rigorous testing and validation of their accuracy. Simulation results demonstrate the precision of the 
developed models, consistently aligning with expected vehicle behavior under different input conditions. Specifically, the models 
accurately replicate vehicle motion in straight lines and circular patterns, corresponding to longitudinal speed and steering angle 
inputs. Additionally, a path tracking controller is integrated to showcase the model's efficiency and validate its derived parameters. 
The reliability and accuracy of the thoroughly developed models underscore their suitability for algorithm development and vali- 
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dation, essential for advancing autonomous vehicle technology and enhancing vehicle safety and performance. 
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1. Introduction 


Since the last few decades, the automotive industry and the market 
structure has changed in an unprecedented way. There were simul- 
taneous growing needs for vehicle safety, environmental protection, 
and intelligent control. As a result, the use of advanced technologies 
such as computer technologies, virtual reality technologies and in- 
telligent algorithms became widespread in the industry. Thanks to 
these advanced technologies, the automotive industry focused on the 
design and testing of new Advanced Driving Assistance System 
(ADAS) functionalities to enhance not only the safety but also the 
performance and comfort of the vehicle and its passengers. The 
ADAS functionalities include lane keeping [1] but also adaptive 
cruise control [2], lane assistance [3], emergency braking [4], and 


parking assistance [5] among others. In recent times, fully Autono- 
mous Driving (AD) is also becoming an increasingly interesting 
topic of research in the automotive industry [6] in which more ad- 
vanced algorithms are required for trajectory planning and control- 
ling. 

Vehicle dynamics plays an important role in the automotive in- 
dustry development [7]. The mathematical models not only allow to 
understand the behavior of the vehicle, but one can also use these 
models to verify and validate the developed algorithms in simulation 
before implementing them on the real hardware [8]. This allows a 
significant cost reduction. Moreover, by studying the complex inter- 
actions between a vehicle's components, such as its tires, suspension, 
and steering system, researchers can identify potential safety hazards 
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and develop innovative solutions to improve vehicle safety [9]. 

Initially, research on vehicle dynamics focused on analyzing the 
vehicle's performance under various external conditions and service 
requirements. During the 1930s, there was a notable shift in research 
focus towards the steering, suspension mechanics, and driving sta- 
bility of the vehicles [10]. Researchers like Lanchester Maurice and 
Segel began to examine how external environmental factors, such as 
road surface roughness, air flow, tire conditions, and driver behavior, 
impacted the vehicle dynamics. Additionally, they studied the inter- 
related effects of these various conditions [11]. In 1993, they pro- 
vided a detailed overview of the progress made in vehicle dynamics 
research prior to 1990 in the Proceedings of Institution of Mechani- 
cal Engineers [12]. 

Over the subsequent decades, significant research has been de- 
voted to investigating vehicle ride comfort and handling stability. 
The study of handling dynamics concerns primarily the lateral or 
transverse dynamics of the vehicle, including its handling, stability, 
vehicle sideslip resulting from tire lateral force, yawing, and roll mo- 
tion. The research on vehicle handling stability in vehicle dynamics 
has evolved from experimental studies to theoretical analysis and 
from open-loop to closed-loop systems. The representative mono- 
graphs of vehicle handling dynamics include "Vehicle Handling Dy- 
namics Theory and Application" by Abe [13] and "Overview on Ve- 
hicle Dynamics" by Yang [7]. 

Vehicle driving dynamics is divided into two categories: longitu- 
dinal and lateral dynamics. Longitudinal dynamics concerns the mo- 
tion of a vehicle along its longitudinal axis, which includes acceler- 
ation, braking, and the resulting changes in velocity. It also considers 
various resistive forces which influence the longitudinal movement 
of the vehicle [14]. Lateral dynamics on the other hand concerns the 
motion of a vehicle perpendicular to its longitudinal axis, which in- 
cludes turning, cornering, and the resulting changes in direction. 
Both longitudinal and lateral dynamics are critical for vehicle han- 
dling, safety, and overall performance. Relevant monographs in this 
area can be found in "Vehicle Dynamics and Control" by Rajamani 
[15]. Another pivotal aspect of vehicle dynamics is the vertical mo- 
tion regulated by the vehicle suspension system. Therefore, a com- 
prehensive vehicle suspension mathematical model is of utmost im- 
portance to analyse and optimize the vehicle’s behaviour. This 
model not only profoundly influences the overall dynamics of the 
vehicle but also serves as a key component of vehicle stability, ride 
comfort, handling, and safety [16]. In recent times, researchers have 
also developed mathematical models not only for electric and hybrid 
powertrains and propulsion systems[17,18] but also for braking 
[14]and speed transmission systems [19]. However, it's crucial to 
highlight that vehicle motion dynamics models serve as pivotal com- 
ponents for testing and verifying these intricate models. 

The main objective of this work is to create precise and efficient 
simulations by delving into key aspects of vehicle behavior. These 
aspects encompass lateral vehicle dynamics, longitudinal vehicle dy- 
namics and vehicle suspension dynamics. The lateral vehicle dy- 
namics pertains to the vehicle’s turning behavior, providing insights 
that enable the effective design of control systems. Longitudinal ve- 
hicle dynamics have been modeled to investigate acceleration, de- 
celeration, and speed control responses to driver inputs and external 
factors. Furthermore, a vehicle suspension systems model has been 
developed to simulate the vehicle’s reaction to uneven road surfaces, 


therefore, enhancing ride comfort and stability, that are integral to 
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Behavior (Lateral, 
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Suspension Dynamics) 


Validate Results of Path 
Tracking Control 
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Environment) 


Figure 1. Flowchart summarizing the entire study. 


passenger safety. Compared to the previous literature, this paper not 
only advances the understanding of the vehicle’s motion but also 
combines all three dynamics models to obtain a better analysis of the 
vehicle’s behavior in real driving scenarios, emphasizing safety con- 
siderations. A detailed derivation of these dynamics is provided in 
the paper. In the next step, these dynamics are implemented in 
MATLAB/Simulink. The implemented models are tested with vari- 
ous input signals. Finally, a path tracking controller is designed and 
implemented. The results of path tracking control show that the de- 
veloped mathematical models mimic the real motion of a vehicle and 
can be used to design and validate algorithms for ADAS and AD in 
a simulation environment. 

The flowchart as shown in Figure | begins with the main objective 
of the study, which is to create precise and efficient simulations by 
investigating key aspects of vehicle behavior (section 1). It outlines 
the key aspects, including lateral, longitudinal, and suspension dy- 
namics (section 2). The detailed dynamics are derived, implemented 
in MATLAB/Simulink, and tested with various input signals. A path 
tracking controller is designed and implemented (section 3), fol- 
lowed by the validation of results in a simulation environment (sec- 
tion 4). 


2. Vehicle model description 


This section discusses the derivation of vehicle mathematical 
modeling. First, the basic coordinate representation is described. 
Then both kinematic and dynamic models are discussed in detail. 


2.1 Coordinate system 


In vehicle dynamics, understanding coordinate systems is crucial 
for accurately analyzing the vehicle’s behavior. In order to develop 
the model, two coordinate frames are required i.e. the inertial frame 
and the body-fixed frame. The inertial frame is a reference frame that 
is fixed relative to the Earth's surface and does not move, while the 
vehicle frame is a reference frame that is fixed relative to the vehicle 
and moves with the vehicle's motion. Figure 2 depicts the reference 
coordinates. The parameters, which are the variable components es- 
sential in deriving the characteristics of the vehicle, are described 
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beneath each corresponding equation. 
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Figure 2. Reference coordination. 


2.2 Kinematic model 


The kinematics of a 4-wheel vehicle can be simplified as 2-wheel 
bicycle model by merging the front and rear wheels. This reduces 
the complexity of dealing with 4 wheels and 2 steering angles to only 
2 wheels and | steering angle. Let us consider that the body reference 
frame is at the center of gravity of the vehicle as shown in Figure 2. 
w describes the orientation of the vehicle, V is the velocity at the 
center of gravity and is the slip angle. 


The motion of the vehicle can then be described as in the Eq. (1) 
below: 


X=V cos(w + B) (1) 
Y=V sine(y +B) 
While the angular speed w can be given as: 
: V 
ee 2 
P= 5 (2) 


Where R is the perpendicular distance from the reference point to 
the Instantaneous Center of Rotation (ICR) represented in Figure 3 
as O. 

Using basic angle geometry, it can be found that the intersecting an- 
gle between the front and rear perpendicular lines is equal to 6, thus 
leading to the following Eq. (3): 


Figure 3. Two-wheel bicycle model [20]. 


ly + Ly i: 
,tan(B) = — (3) 


TR TR 


tan(d) = 


Where J, and ly are the distances between the center of gravity 
and the rear and front axles respectively, rg represent the distance 
between O to center of rear axle. From the above equations, the fol- 
lowing Eq. (4) is obtained: 


wp = +5 cos(B) tan(6) 


: (4) 
= t a 
B = arctg ( Dae 


tan(6)) 


Eq. (1) and Eq. (4) combined describe the kinematics of the vehicle. 


2.3 Dynamic model 


In order to keep the model simple, some key assumptions were 
made. First, the vehicle is considered as a rigid body, whose dynam- 
ics is determined by the fundamental laws of motion. Second, the 
steering angle is small, which allows us to make the small-angle ap- 
proximations such as cos(d) ~ 1 and sin(d) = 0. 


" Lateral Dynamics 
Let us consider the free body diagram shown in Figure 4. Accord- 
ing to Newton’s second law: 
YF=ma and YM=Iz (5) 


Where }'F represent the sum of forces, m represents the vehicle 
mass and a represent the vehicle acceleration. As well as )} M_ rep- 
resent the sum of all moments, Iz represents the moment of inertia 
and ¥ represents the angular acceleration of the vehicle. 

Considering the lateral forces on front and rear tires in the small- 
angle approximation. The lateral acceleration is a combination of V,, 
which is due to the motion along the y axis, and YV,, referred to as 
centripetal acceleration. Therefore, the following Eq. (6) is obtained: 

Fyp + Fy = m(V, + ¥Y,) (6) 


Now taking the moment around the center of gravity G: 
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Y 


Figure 4. Free body diagram [20]. 


Fyply — Fyrl, = 18 (7) 
The slip angle of front and rear tire can be expressed as: 
YWU+l@ VY, -L@ 
¢=6= and eg, = 2 (8) 
Vr ro Vy 
Using the cornering stiffness of the tire, F,- and F,, are expressed 


as in the following Eq. (9): 


Fy = Cy Op and FE 


y ‘yr = —C,a, 


(9) 


Combining all above equations, the lateral model in the state- 
space can be represented as in Eq. (10): 


‘ 0 1 0 0 
Crtcr Cel ¢—Crly 
le) | ° ae Ot 
a|p| — | 0 0 0 1 
j _ eflp-erlr - cp +lp cy 
¥ 0 IZVx e IzVx 
0 
vy] {2 
x Ly| + fg (10) 
i 
vg I 


" Longitudinal Dynamics 


The longitudinal dynamics describes the motion and behavior of 
vehicles along their longitudinal axis, which involves acceleration, 
deceleration, braking, and the associated forces and interactions af- 
fecting the vehicle's motion. It directly influences the overall perfor- 
mance and handling characteristics of a vehicle. The ability to accel- 
erate smoothly, maintain stability during braking, and ensure precise 
control over the vehicle's longitudinal motion are essential for acom- 
fortable and safe driving experience. Consider a vehicle moving on 
an inclined road, the longitudinal forces acting on the body are 


10 


shown in Figure 5. 


Figure 5. Longitudinal forces acting on the vehicle [15]. 


While in motion, a vehicle encounters multiple resistive forces 
due to various factors acting in the opposite direction to the move- 
ment. These resistance forces include rolling resistance, aerody- 
namic resistance, and weight resistance [14] . The force balance 
along the longitudinal axis leads to the following Eq. (11): 
mx= Fyxft+F xr—Faero—Rxf¢—Rxr—mg sind (I 1) 
Where F'xf represent the longitudinal tire force at the front tires, Fxr 
represent is the longitudinal tire force at the rear tires, Faero repre- 
sent the longitudinal aerodynamic drag force, Rxf represent the the 
force due to rolling resistance at the front tires, Rxr represent the 
force due to rolling resistance at the rear tires, m represent the mass 
of the vehicle, g represent the is the acceleration due to gravity and 
6 represent the angle of inclination of the road on which the vehicle 
is traveling. 

The equivalent aerodynamic drag force on a vehicle can be repre- 
sented as [21]: 


1 
Fyero = 5 PCaAr Ve + Vigina (12) 


Where, p is the mass density of air, Cg is the aerodynamic drag co- 
efficient, A; is the frontal area of the vehicle, V, is the longitudinal 
vehicle velocity and Vying is the wind velocity. 

Similarly, the tire force can be expressed as in the following Eq. (13): 
Pyp = Caf Axp and Fey = Cay xr (13) 
Where Cy¢ and Cg, represents the longitudinal tire stiffness pa- 

rameters of the front 


and rear tires respectively, and @, represents the slip or tire skid- 
ding and given as: 


Tref Ow — Vy 


Ay = 
Vy 


(14) 


Where R,, is the tire rolling resistance which refers to the force that 
opposes the motion of the tires as they roll on the road surface. If the 
road is inclined at an angle a, the rolling resistance R, can be ex- 
pressed as: 


R, = C,,mgcos(a) (15) 
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Where C,,, is the rolling resistance coefficient. It depends on the tire 
construction, the materials, the air pressure, the vehicle speed, and 
the road conditions. In addition, m is the vehicle mas and g rep- 
resent the gravitational acceleration. 


= Suspension Dynamics 


The suspension system plays an important role in the stability of 
a vehicle operation [22]. It stands as a pivotal element that pro- 
foundly influences the vehicle dynamics and passenger comfort. 
Suspension systems can be categorized in three main classes: passive, 
semi-active and active. Passive suspension systems are traditional 
and common suspension systems that mainly rely on mechanical 
components like coil spring, shock absorber and connecting ele- 
ments [16]. On the other hand, active suspension systems use elec- 
tronic sensors and actuators to actively control vehicles height and 
damping characteristics. Semi-active suspension systems bridge be- 
tween passive and active type by using just electronically control 
dampers. 
In this paper, the work is focused on passive suspension systems. 
Figure 6 shows a schematic of a 2 Degree of Freedom (DOF) quarter 
car suspension. 


Vehicle body’ 


ca 


22 


Figure 6. Two-degree of Freedom (DOF) quarter car suspension. 


The equations of motion related to mass m, and my are given in 
Eq. (16) and Eq. (17) by: 


_ : ; 16 
m2, = —k,y(z — 22) — (4, — 22) — Fy (16) 


Moz, = ky (4 — 22) + (4 — 42) —kn(z, — u(t) — Fh (17) 
Where m, is the body mass, mz is the wheel mass, k,& kzare 
the stiffness coefficients, c,is the damping coefficient, z, is the dis- 
placement of the body mass, Z zis the displacement of the wheel 
mass, F, is the applied external force and u(t) is the road input 
displacement. The state-space can be obtained from Eq. (16) and Eq. 
(17) such as: 
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Zsf, 


A. 
Unspring mass 


=a 


: 0 1 0 0 
% (2 a aS 2 
Z| _ m4, my, my, my, 21 
Z| 0 0 0 g- 1? la” 
Zo Ky (ka +k) C1 Zo 
M2 m2 m2 m2 (18) 
0 0 
A. a 
m4 24 
0 0 | bal and y(j= [1 0 1 O] Zy 
Fe Zp 
m2 ™g2 


A more comprehensive model can be represented through a half car 
suspension model. It considers a half car, including front and rear 
body parts. The free body diagram represented in Figure 7 is consid- 
ered. 
With the small angle approximation, the equation of motion re- 
lated to half body mass m, is given from Eq. (19) by: 
M,Z, = —k (zs =f ad) —k,(Z5 — Zyr + bO) 


— ¢,(Z; — Zur + DO) 


Where: m, is the body mass, @ represents pitch angle, ky, k,, 
Kz, Ker are the stiffness coefficients, a and b are the distances 
from the front and rear axle to the center of gravity respectively, Zz, 
is the displacement of the body mass, Z,¢ and Z,, are the displace- 
ments of the front and rear body masses respectively, Z,¢ and Zy; 
are the displacements of the front and rear wheel masses respectively. 

While the equations of motion related to front and rear wheels, 
Myf and m,; respectively, are given from Kq. (20) and Eq.(21) as: 


Muphup = Kp (Zs — Zug — AO) — ker (Zu — Zr) 
(20) 


+ cr (2 - Zur _ ad) 
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MurZur = k, (Zs — Zyyp + bQ) 7 Ker (Zur ~ Zrr) 
(21) 

+ Cy(Z, — Zur + DO) 
Where m,, and m,, are the front wheel and rear wheel masses 
while cr and c, are the damping coefficients. While z,- and 2, 
are the road input displacements for the front and the rear wheel re- 
spectively. 

In the half car model, the moment around the center of gravity 
must also be considered, such as in (22): 


1,0 = ak, (2, — Zug — AD) — Dk (Z, — Zur + BB) 


+ acs (Zs — Zur — aD) (22) 
— be,(z5 — Zyr + bD) 
Where Jy represents the moment of inertia. 
The following state matrix is then obtained as in (23): 
=| ae O ny Ge ZO 2y Ze J 3) 


The state space represented here after is then obtained as in (24) — (28): 


000010 0 0 0 0 
000001 0 0 0 0 
x=|0 0 0 0 0 0 1 O}x+]0 Of} u (24) 
0000000 1 0 0 
—-M1K —-M"C M1F 
Where: 
Cp +c; bc, — acp Cf Cy 
Ce be, —acg cra*+c,-b? ace —be, (25) 
=o aC Cr 0 
—C, —bc, 0 Cr 
m, 0 O 0 
0 lL oO 0 26 
Me 0 O my 0 9) 
0 0 O m, 
K= 
kp +k, bk, — akg —Ky —k, 
bk, —akp ka? +k,-b?  akp — —bk, @)) 
—K ak Kp + ker 0 
—k,. —bk, 0 k,. + Ker 
0 0 
hae (28) 
~ [ker 0 
0 Key 
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Figure 8. Representation of the principle of the Pure-Pursuit control. 
3. Path tracking 


A Pure-Pursuit control algorithm for lateral control was imple- 
mented with the developed dynamics of the vehicle in order to vali- 
date the model. In this section, a brief explanation of this controller 
is presented. 

The pure pursuit technique [23,24] involves a geometric compu- 
tation of the curvature of a circular arc, linking the position of the 
rear axle to a target point situated along the path ahead of the vehicle. 
This objective point is established based on a look-ahead distance 
lq extending from the present rear axle position towards the in- 
tended path. The target point (g,, gy)is illustrated in Figure 8. 

The steering angle 6 of the vehicle can be found using solely the 
coordinates of the target point and the angle o formed between the 
vehicle’s heading vector and the look-ahead vector, leading to the 
following equation: 


(29) 


( 2L —_) 
6 = arctan| —————— 


la 


4. Simulation results and discussion 


The vehicle models for both longitudinal and lateral dynamic are 
implemented in MATLAB/ Simulink as shown in Figure 9. Longi- 
tudinal model of the vehicle computed through (11)-(15) is imple- 
mented in ‘longitudinal vehicle dynamics’ block while lateral model 
represented by (10) is realized in ‘lateral vehicle dynamics’ block. 
For the simulation purpose, the parameters of a Tesla Model S were 
considered. By incorporating these specific parameters, the objective 
is to conduct comprehensive simulations that closely mimic the real- 
world scenarios. 

A half-car suspension model is also’ simulated in 
MATLAB/Simulink. The results are illustrated in Figure 10 (a-b). 
The suspension system is subjected to two different disturbance sig- 
nals for the front and rear wheels. These results demonstrate the evo- 
lution of body displacement. An examination of the graphs reveals a 
key principle of this system’s operation: the maximum body dis- 
placement is half that of the wheel’s maximum displacement, high- 
lighting an essential characteristic of the system’s behavior. 


F. Haidar et al. 


Engineering Perspective 4 (1): 7-16, 2024 


Vehicle longitudinal and lateral dynamics are tested with various 
inputs. The initial states of the vehicle are set as Xp = 0, Yo = 0, 
Wo = pi/4 and V,, = 0m/s. First, the steering angle is fixed as 0 
with a lateral speed of 20km/h. A Proportional-Integral-Derivative 
(PID) controller is used to achieve desired lateral speed. The tunning 
proportion, integral and derivative tunning gains are chosen as 20, 4 
and 3 respectively. Figure 11 shows the obtained results for this case. 
Figure 11-a depicts the longitudinal speed while Figure 11-b illus- 
trates the position of the vehicle. Since the steering angle is fixed, 
the vehicle moved in a straight line. In the second case, the steering 
angle of the vehicle is kept constant at the value of 1/16 rad 
while keeping the same longitudinal speed profile and initial condi- 
tions. The obtained result is shown in Figure 12. As anticipated, the 


vehicle executed a precise circular motion, showcasing the accuracy 
of the model. 

The path tracking task involved following a predefined path repre- 
sented by a series of waypoints. The pure-pursuit controller is used 
to control the lateral motion of the vehicle through appropriate steer- 
ing commands. The lookahead distance was dynamically adjusted to 
ensure a smooth tracking. Figure 13 shows the path tracking results 
along with the trajectory error. One can see that the maximum track- 
ing error is always less than 0.025 meter which indicates that the d 
esired trajectory is followed quite accurately. Moreover, the simula- 
tion results also signify the accuracy of the developed mathematical 
models and controller. 
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Figure 9. Detailed Simulink model for longitudinal and lateral dynamics. 
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Figure 11. Longitudinal speed (a) and position of the vehicle (b) with a steering angle of zero. 


5. Conclusions 


In this paper, a mathematical model of vehicle dynamics is devel- 
oped to describe a vehicle’s motion in lateral, longitudinal and ver- 
tical directions. These dynamics models are vital for developing and 
validating algorithms for the vehicles. The models are implemented 
in MATLAB/Simulink and simulation are carried out with various 
input signals. The simulation results clearly validate the precision of 
the developed models, as they consistently mirror the expected ve- 
hicle motion. For instance, execution of motion by the vehicle in a 
straight line or in a circular pattern according to the longitudinal 
speed and steering angle inputs, the simulation outcomes align seam- 
lessly with the anticipated behaviour. These results underscore the 
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reliability and accuracy of the meticulously developed models, af- 
firming their efficacy in capturing and replicating the dynamic re- 
sponses of the vehicle in various scenarios. Moreover, a path track- 
ing controller is also integrated to demonstrate the efficiency of the 
model and verify the accuracy of its derived parameters. Integrating 
a powertrain model with the vehicle dynamics is anticipated as a fu- 
ture work. 


F. Haidar et al. 


Engineering Perspective 4 (1): 7-16, 2024 


25;-——— 


Error Plot 


Path and Trajectory Plot 


20 (a) 


= 
a 


Y-position 
3 


Initial position 


a 


-10 5 0 5 
X-position 


-20 -15 


Figure 12. Position of the vehicle with a steering angle of 7/16. 
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